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1. Introduction 

In a recent paper, Dennis, Kitaev, Landall and Preskill (DKLP) 1 showed that error cor- 
rection in the toric code for stable storage of quantum information is closely related to the 
phase transition between ferromagnetic and paramagnetic phases in the two-dimensional =b J 
Ising model of spin glasses. This is remarkable in that it paves a way for transferring results, 
concepts and insights between the two at first sight unrelated domains of quantum error 
correction and spin glasses, and provides the motivation for the present work. 

An important aspect of the toric code is that it has a finite error threshold; this is the 
critical value of error probability per qubit beyond which it is impossible to correct errors in 
the thermodynamic (large system-size) limit. Thus, for error rates larger than the threshold, 
the encoded quantum information is lost. It is therefore necessary to estimate the precise 
value of the error threshold for, e.g., the design of hardware of quantum memory. DKLP have 
shown that this error threshold is equal to the probability of antiferromagnetic bonds at the 
multicritical point in the phase diagram of the ± J Ising model on the square lattice. 

We do not discuss the value of error threshold itself in the present paper. 1-3 We instead 
strive to clarify a closely related problem of the structure of error patterns underlying a 
characteristic feature of the toric code, degeneracy, by deriving bounds on the number of 
equivalent error patterns. These bounds are shown to be related to the information-theoretical 
entropy of the distribution of frustrated plaquettes in the ± J Ising model. It will also be shown 
that a lower bound on the ground-state energy of the ± J Ising model can be obtained from 
the bounds derived for the toric code. 
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We outline the link between spin glasses and quantum error correcting codes in §2. Errors 
are detected by measurement of a syndrome, from which the underlying error pattern is to be 
inferred and then corrected. Due to the nature of quantum encoding in the toric code, however, 
many error patterns are equivalent (degenerate), and it is sufficient to infer the equivalence 
class of the true error pattern, not the precise pattern itself. In §3 we count the number of 
syndromes D and the total number of equivalence classes C, which is an easy task. The hard 
part is counting the number of equivalence classes C(p) containing error patterns with a given 
fraction p of errors. We first ignore the issue of degeneracy, and count error patterns instead 
of equivalence classes, from which upper bounds on C(p) are derived in §4.1. Lower bounds 
are discussed in §§4.2 and 4.3. These bounds are summarized in §4.4. One of the lower bounds 
involves the ground-state energy of the =L J Ising model, and this fact is used in §4.5 to derive 
a lower bound on the latter ground-state energy. The results are summarized and discussed 
in §5. 

2. Toric code and spin glass 

It is useful to first sketch the connection between toric quantum codes and two-dimensional 
spin glasses, following DKLP, 1 since most of the readers may be unfamiliar with this relatively 
new interdisciplinary field. 

2.1 Encoding in the toric code 

Consider a square lattice with N = L 2 sites and n = 2N bonds labelled by b = (ij) 
and connecting nearest-neighbour sites i and j; toroidal boundary conditions are used so that 
bonds on a boundary connect to the opposite side of the lattice. The toric code comprises n 
quantum (spin-1/2) spins located on the bonds of the lattice. We call the local Pauli operators 
Xf,, Yb and Zb- In the Z-basis of the state space, each basis vector \z) = \{zb}) is specified by the 
z-components Zb = ±1 of all spins. The Pauli operators then act as Z\\ . . . Zb . . .) = zj,\ . . . zj, . . .) 
and Xb\ . . . Zb ■ ■ ■) = \ ■ ■ ■ —Zb ■ ■ ■) so that Xb effectively just flips the spin at b; also, Yj, = iXbZb- 

Quantum states are vulnerable to decoherence, and robustness should be introduced to 
protect quantum information. Quantum error correction is a powerful method for this purpose, 
in which one encodes quantum information by mapping it to another (redundant) set of 
quantum states. Specifically, in the toric code, one maps the state space of two logical qubits 
that are to be encoded onto a 2 2 = 4-dimensional subspace (the "code space") of the spin 
system's 2 n -dimensional state space. This mapping is redundant because n qubits are used to 
represent 2 qubits. A basis \tpo), IV^) for the code space can be defined as 



zed 

where Cj denotes a class of states (equivalence class) to be defined shortly. The sum here runs 
over states \z) which form cycles within an equivalence class. This means that out of the four 




(1) 
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Fig. 1. Example of a cycle. Full lines represent bonds with Zb = 1 and wavy lines are for zi, = — 1. 
Dual bonds to z b = — 1 form a closed loop (shown dotted). 



bonds b(P) around each lattice plaquette P an even number are negative (zb = — 1). A cycle 
\z) is thus an eigenstate with eigenvalue 1 of all the operators Zp = Y\b(p) Zb- Geometrically, 
\z) is a cycle if the duals to its negative bonds form closed loops on the dual lattice (see Fig. 1). 

2.2 Equivalence class 

Cycles are called equivalent if they can be locally deformed into each other, by repeatedly 
flipping all spins around some plaquette of the dual lattice. As illustrated in Fig. 2, this 
corresponds on the original lattice to applying a product of operators of the form Xj = 
Y\b(j) ^bi where the b(j) are the four bonds meeting at site j. It is then easy to see that there 
are exactly four equivalence classes of cycles, denoted Co,---,Cz: Co contains all the trivial 
cycles, which are equivalent to the empty cycle state \{zb = 1}). C\ and C2 collect cycles 
equivalent to a single loop winding across the lattice boundary in the horizontal and vertical 
directions respectively, and C3 those with both a horizontal and a vertical loop (Fig. 3). In 
(1), what is meant is that for each |V>i) the sum runs over all cycle states \z) in the equivalence 
class C{. This implies in particular that the are invariant under the action of any of the 
operators Zp and Xj: the Zp leave each cycle state \z) invariant, and the Xj only permute 
cycle states within an equivalence class, 

Z P \^ l )=X^ l ) = \^ l ). (2) 

In the large-A limit, the toric code has zero code rate R: it encodes k = 2 qubits using 
n = 2N qubits, so that R = k/n = 1/N — > 0. The point of this highly redundant encoding is 
to allow for the correction of errors arising from decoherence caused by the interaction of the 
quantum state with its environment. An error introduced in this way corresponds to a product 
of Pauli operators X, Y, Z acting on some of the quantum spins. Because of Y = iXZ, Y- 
errors can be treated as combinations of X- and Z-errors. We write a state with A-errors on 
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Fig. 2. Four equivalent cycles. Full, wavy and dotted lines have the same meaning as in Fig. 1. 
Application of X± to (a) gives (b); similarly, (c) and (d) are obtained from (a) by applying X 2 Xi 
and X3X2X1, respectively. 



a set of bonds S x and Z-errors on S z as 

m = ( n x >> n ^ ) i^>- (3) 

If the product of and Z& on the right-hand side can be represented by a product of Xj and 
Zp only, then = \ipi). In general, however, this is not the case, \ipi) ^ \ipi). 

Furthermore, the toric code is in the general class of Calderbank-Shor-Steane (CSS) 
codes, 4 ' 5 for which the X- and Z-errors can be treated separately, without interference be- 
tween the corresponding error correction procedures because [Zp,Xj] = for any P and j. 1 
We can therefore focus in the following exclusively on X-errors, i.e. spin flips. Z-errors can be 
discussed separately in the same manner but on the dual lattice. 1 

2.3 Syndrome and error correction 

We define an error pattern {r^} such that 17, = — 1 if an X-error has occurred on bond b 
(i.e. if b 6 S x ) and 17, = 1 otherwise. To diagnose where errors have occurred, one measures all 
the Zp. Without corruption all measured values would be 1 according to (2). The errors give a 
nontrivial set of measurement values Hfe(p) Tfe = — 1 at those plaquettes around which an odd 
number of errors have occurred, giving the syndrome, the set of plaquettes with measurement 
value —1. It should be stressed that the syndrome measurement does not cause any quantum 
decoherence: From (3), one easily sees - by commuting Zp to the right through the various 
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(a) 



(b) 




(c) 



(d) 



Fig. 3. Representative cycles of four different classes, Cq,C\,C?, and C3 in (a), (b), (c) and (d), re- 
spectively. 
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Fig. 4. Error patterns with the same syndrome (crosses). Double lines are bonds where Tb 
Error chains are written in dashed lines. 



Xf, and using (2) - that even a corrupted state is an eigenstate of each Zp. Explicitly, 
Zp\ipi) = with the plus sign when there are an even number of errors around P and 

the minus sign otherwise. 

One can visualize the error pattern r by drawing the duals to the bonds with T& = — 1; 
these form "error chains" ending in "defects", i.e. plaquettes where the syndrome has detected 
an error (Fig. 4). The syndrome measurement is therefore highly ambiguous: any error chain 
t' with the same set of defects as r gives the same syndrome as exemplified in Fig. 4. The 
condition for this is n&(p) T b = rib(p) T 'b f° r au plaquettes P and hence n&(p) T b T b = 1 ; the 
bonds with t\t[ = — 1 form a cycle rr' in the sense defined above. 
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Fig. 5. Non-equivalent error patterns with the same syndrome. 



Now assume we have inferred some error pattern r' consistent with the syndrome (e.g. 
Fig.4 (a)), so that rr' is a cycle, where r is the actual error pattern (e.g. Fig.4 (b)). If we 
correct errors according to r', by applying a spin-flip to all spins with r' b = —1, this 
can be viewed as first correcting the errors u that actually occurred, followed by a series of 
spin-flips where r^ = — 1 (four double lines in Fig. 4 (a) and (b)). The first stage recovers 
the uncorrupted code state \if)i). If tt' is a trivial cycle (rr' € Co), then the second stage 
corresponds to applying a product of operators of the form Xj. But these leave code states 
invariant, see (2), so that our error correction was successful. This is a key difference of this 
degenerate quantum code to classical error correction: we do not need to detect all details of 
the error pattern r, but only its equivalence class (defined as the set of all r' such that rr' 
is a trivial cycle). Error correction according to r' will be unsuccessful, on the other hand, if 
r and r' are non-equivalent, i.e. if rr' is a nontrivial cycle like Fig. 5: the second stage from 
above then mixes up the basis vectors of different code spaces. 

2.4 Spin-glass representation 

We can now relate the decoding task to a spin glass problem. The condition that rr' is a 
trivial cycle can be written as r|- = TijOiOj if we revert to denoting each bond b by its end 
points i and j. Here the <7j = ±1 are classical spin variables. Then the second (conceptual) 
stage of error correction discussed above would consist of the application of an operator Xj 
for each site where Oj = —1 (the central site in the case of Fig. 4). Assuming that errors 
were generated independently and with probability p on each bond b, we assign an overall 
probability 

p{j') = n(i -pf+^/y 1 -^/ 2 = b(i -P)r /2 ex P UpXH ) ( 4 ) 

(«) \ (ii) / 

to any error chain r'; here we have defined K p = ^ln[(l —p)/p]. The total probability for 
error chains r' which lead to successful error correction is therefore proportional to 

Z = ex P \ K P T>i A = Y exp K P T ij a i a J ■ ( 5 ) 

t'-.tt'gCo \ (ij) I a \ (ij) / 
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This is the partition function of an Ising ± J spin glass with local interaction strengths Jij = 
KpTij] the relation between K p and p implies that the system is on the so-called Nishimori 
line (NL). 6 ' 7 If the actual error pattern r was indeed generated according to the assumed 
probability weight (4), each bond in (5) is ferromagnetic with probability 1—p (jij = 1) and 
antiferromagnetic otherwise. 

The total probability for error chains r' which lead to faulty error correction has the same 
form as (5), but with modified boundary conditions. For example, any r' such that rr' G C2 
can be written as t-j = TijOiOj with the convention that all bonds along, say, the left boundary 
of the lattice are inverted. This corresponds to the use of antiperiodic boundary conditions in 
the left-right direction when evaluating the spin products OiOj. The total probability weight 
Z<i for error chains t' with rr' G C2 is therefore a partition function of the form (5) with these 
modified boundary conditions; in Zq we had implicitly assumed periodic boundary conditions. 
Similarly Z\ and Z% correspond to the partition functions for antiperiodic boundary conditions 
in the top-bottom direction, and in both left-right and top-bottom directions. 

Combining the above results, we conclude that our total probability of inferring from the 
syndrome an error chain r' which leads to successful error correction is Zq/(Zq + Z\ +Z2 + Z3) 
and close to unity as long as Zq 3> Z^ for k = 1,2, 3. This condition is met if p is small enough 
so that we are in the ferromagnetic phase of the spin system defined by (5): the existence of 
domain boundaries then implies for the free energies Fi = —TlnZi that Fk — Fq {k = 1, 2, 3) 
is positive and of order L, thus Zq 3> Z\, 2*2,2%. In the paramagnetic phase, on the other 
hand, — Fq = Oil) and there will be a nonzero probability for error correction to fail. 
In summary, the toric code can correct errors below an error threshold, i.e., in the range 
< p < p c where the associated ± J Ising spin glass on the NL is in its ferromagnetic phase 
on the square lattice; p c is then the location of the multicritical point. We may therefore be 
able to learn something about the spin glass problem from knowledge about the toric code 
and vice versa. This is our motivation for the present work. 

3. Simple number counting 

The argument in the previous section suggests that the numbers of possible syndromes 
and equivalence classes of error patterns would give important measures of performance of 
error correction in the toric code. We therefore discuss this problem in the present and next 
sections. 

Let us first count the total numbers of different syndromes and equivalence classes, without 
specifying the error probability p. Since at each plaquette the syndrome measurement of Zp 
can give either +1 or —1, the total number of syndromes is 

D = 2 N ~\ (6) 

using that for the square lattice the number of plaquettes is equal to the number N = I? 



7/15 



J. Phys. Soc. Jpn. Full Paper 

of lattice sites. The factor 2~ 1 arises because only an even number of sites with nontrivial 
syndrome Zp = — 1 can exist, as can be seen from the fact that the product of all Zp is the 
identity operator on the torus (all Zf, appear twice from neighbouring plaquettes and Z? = 1). 

The total number of error patterns is 2™. And the number of error patterns in an equiv- 
alence class is 2 N ~ 1 because, as explained above, equivalent error patterns are related by 
r-j = TijGiGj. Each of the Oi = ±1 can be chosen independently and gives a different r', 
except for an overall reversal of the spin configuration which leaves r' unchanged and gives 
the factor 2 _1 . Thus the number of equivalence classes is 

C = 2^1 = 2 " +1 ( 7 ) 
which is equal to 4D. This correctly indicates that each syndrome corresponds to four different 
equivalence classes of error patterns. The above argument easily generalizes to other lattices: 
in general, D = 2 P ~ 1 where V is the number of plaquettes, while C = 2 n /2 N ~ 1 = 2 n - N+1 . 
The equality C = AD then follows from Euler's theorem n = V + N . 

If we now specify the fraction of errors or error probability p, most of the 2 n error patterns 
are excluded because the number of errors these patterns have is different from np. Thus the 
number of equivalence classes containing patterns with np errors, C(p), is significantly smaller 
than C. Then the number of syndromes D (= 2 iV ~ 1 ) is much larger than C(p), and we have 
a sufficient number of syndromes to specify the equivalence class corresponding to a given 
syndrome, D = C/4 S> C(p). This implies that a simple number counting does not lead to 
the critical value p c by the classical argument that, beyond p c , the number of errors exceeds 
that of syndromes (C(p) > D) and one cannot identify the error from the syndrome, leading 
to unsuccessful error correction. 

4. Bounds on equivalence classes and ground-state energy 

We next proceed to evaluate C(p). Upper and lower bounds are derived for this quantity, 
which will further be shown to give a lower bound on the ground-state energy of the ± J Ising 
model. 



4-1 Upper bounds 

Let us begin by constructing upper bounds. A trivial upper bound for C(p) is C, 

C{p) <C = 2 N+l w 2 N , (8) 

where the last expression gives the result to exponential accuracy, which is all we are normally 
interested in. To derive another upper bound, let us temporarily ignore the equivalence of error 
patterns, which leads to overcounting. The number of classically distinct error patterns with 
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a fraction p of errors is 

~2 nH( - p \ (9) 




where H(p) is the binary entropy, 

H(p) = -p log 2 p-(l-p) log 2 (l -p). (10) 

We therefore have 

C(p) < 2 nH ^ (11) 

because equivalence will reduce the number by grouping errors into classes. The two upper 
bounds (8) and (11) cross each other at the point where H(p) = \. 

An interesting observation is obtained if we continue to ignore the issue of equivalence 
and apply the argument for classical codes to the present problem. The number of syndromes 
is D = 2 Ar_1 rs 2 N . Thus, if we demand that there exist a sufficient number of syndromes to 
distinguish all the errors, we have 

2 nH(p) < 2 JV _ ( 12 ) 

On the square lattice, n = 2N, and (12) leads to H(p) < \ or p < 0.1100 as a necessary 
condition for classical non-degenerate error correction to be successful. This boundary value 
H(p c ) = \ happens to coincide with the conjecture on the exact location of the multicritical 
point (separating the ferromagnetic and paramagnetic phases on the NL) of the ±J Ising 
model on the square lattice, which agrees well with numerical estimates (see 2 ' 3 and references 
therein). Somewhat surprisingly, therefore, a naive argument which ignores degeneracy nev- 
ertheless seems to give the correct value of the error threshold for the highly degenerate toric 
code. 1 

It is also interesting that the lower bound on the existence of generic error-correctable 
CSS codes, if applied to the toric code, coincides with the above result as pointed out by 
DKLP: 1 It is known that there exist CSS codes with critical error probability p c and code 
rate R = 1 — 2H{p c ) in the asymptotic limit of large code size. Since the toric code has R — ► 
asymptotically, one finds H{p c ) = ^. 



4-2 Lower bound (I) 

Lower bounds are derived by slightly more elaborate arguments. A loose lower bound is 

2 nH (p) 

C(p) > (13) 

The denominator on the right-hand side is the maximum number of error patterns in an 
equivalence class, thus leading to a smaller value than C(p). 
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A stronger lower bound is given as 

OnH(p ) 

C(p) > (14) 

where po is a function of p denned by 

E 9 (po) = -n(l - 2p) (15) 

for sufficiently large N. Here E g (po) is the ground-state energy of the ±J Ising model with 
a fraction po of negative bonds in a typical configuration, i.e. one where the positive and 
negative bonds are randomly distributed. To prove (14), we first recall from above that the 
number of error patterns for given p, 



Mp) = E S (E^- E p\ = ( «2^W, (16) 

r J \np J 

where E p = n(l — 2p), is an overcount of the number of equivalence classes because it ignores 
equivalence. Each term in (16) should be divided by the number of error patterns equivalent 
to r which contain the same number of errors: 

B & t ) = \ E 6 E - E v\- (17) 

The factor \ comes from overall up-down symmetry (see before (7)). We therefore have 



.HY,(ij)Tij -E p 

B( P ,T 



The prime in the sum indicates that terms with B(p, r) = should be excluded. Applying a 
gauge transformation r%j — > TijOiOj and averaging over the gauge variables \oi\ gives, up to 
an unimportant factor of ^, 

Bfc) " 2^ E = 237 E °( S (P. ( 19 ) 

where 0(x) is 1 for x > and for x = 0. This result is intuitively clear: the final sum in (19) 
counts all error patterns r for which an equivalent pattern with a fraction p of errors exists. 
Every equivalence class contained in C(p) thus contributes 2^ v ~ 1 ~ 2 N times to the sum, and 
the prefactor compensates for this. 

The lower bound (14) is now obtained by restricting the sum in (19) to typical error 
patterns r with npo errors (i.e. npo negative r^'s) 

°(P) * W E e(B(p,r)). (20) 

r(po) 

For a typical configuration r with npo negative bonds, B(p, r) is almost always positive because 
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the constraint on the right-hand side of (17) 

^2 T ij a i a j ~ E p = (21) 
(ij) 

is satisfied by a ground-state configuration of the Uj, due to the definition of p in (15). We 
therefore find 

C(p)>±( n \ (22) 



npo 



2 N 

which is the lower bound in (14). 
4-3 Lower bound (II) 

Another lower bound is obtained by restricting the sum in (18) to typical bond configu- 
rations r with np negative bonds (errors). Then B(p, r) is - in spin glass terms - the number 
of spin configurations with energy — E p = — n(l — 2p) for a typical configuration r with np 
negative bonds, the logarithm of which is the thermodynamic entropy on the NL: 6,7 

B(p,r) = e 5(p ' T) . (23) 

Since r is a typical configuration, S(p, r) does not depend on the details of r for sufficiently 
large system size, so we denote it as S(p). Thus the denominator on the right-hand side of 
(18) can be brought in front of the sum and we find 

C(p) > e- s ^^5 K>, - E p \ « 2"^)- s W/ In2 . (24) 

t(p) \ (ij) ) 

4-4 Behaviour of bounds 

Our upper and lower bounds can be summarized as 

max (V"(p)-S(p)/m 2) 2 nH( po )-N^ < c{p) < min (^iv 2 nH( P )^ (25) 

This result is depicted in Fig. 6, where we are using for S(p) an upper bound from our previous 
work. 8 It is seen that C(p) reaches its maximum value 2 N at some p below 0.15 and above 
0.1100. This is because one of the lower bounds 2 nH ( p °)~ N reaches the upper bound 2 N at 
p close to 0.15. The latter value was obtained by the relation nH(po) — N = N (or po = ^) 
and using the numerical value of the ground-state energy E g (^) = —1.407V (see Fig. 7 for this 
value) . 

This saturation of the upper bound at an intermediate value of p is not unnatural because 
C(p) is the number of classes of equivalent error patterns, with all classes counted with uniform 
weight. If we instead give appropriate probability weights to various error patterns, and thence 
to equivalence classes, we would reach a smaller value, the logarithm of which we denote by S n . 
Apart from a trivial factor of In 2, this quantity is nothing but the lower bound 2 nH ^~ s ^ p ^ ln2 



11/15 



J. Phys. Soc. Jpn. Full Paper 




P 



Fig. 6. Upper and lower bounds of C(p) from (25), plotted as bounds on TV -1 log 2 C(p). Solid lines 
give the upper bounds, 1 and 2H(p), and dashed lines the lower ones, 2H(p) — S(p)/(N\n2) and 
2H(po) — 1. Because S(p) is not known exactly, we replace it by an upper bound S+(p) from our 
previous study. 8 To sketch 2H(p ) — 1, an rough fit was made to the numerical data in Fig. 7 to 
obtain p(po) and from there the inverse relation po(p)- 



because the latter was derived by using typical error configurations, 

S w = nH(p)]n2-S(p). (26) 

This S n is of course smaller than the maximum value of In C (p) and reaches its maximum 
only at p = \. 

It is instructive to derive (26) from a different argument. Consider the probability weight of 
an error pattern r as given in (4) . The probability weight of the equivalence class it containing 
this error pattern is then 

P K = \p(l-p)] n/2 Zo(p,T), (27) 

where Zq(p,t) is the partition function (5) for the given p (on the NL) and r. Then the 
information-theoretical entropy of the probability distribution P w of equivalence classes is 

S 7r = -J]P 7r lnP 7r , (28) 

7T 

which is equal to the information-theoretical entropy of distribution of frustrated plaquettes in 
the ± J Ising model. 8 ' 9 Using (27) and the relation F = —Tin Z = E—TS with E = —n(l—2p) 
on the NL, it is possible to reduce this expression to (26). To see this, (28) is first rewritten 
using (27) as 

n E 
S7T = --ln[p(l-p)] + --S(p). (29) 
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Fig. 7. Ground-state energy of the two-dimensional ± J spin glass. Shown is the lower bound resulting 
from (30), compared with numerical results for linear system sizes L = 40, 50 and 60. The latter 
were obtained by averaging over 10, 15 and 5 numerical realizations of the disorder, respectively; 
error bars are smaller than the symbol sizes. All numerical calculations were run on the spin glass 
ground state server. 10 



With E/T = -n(l -2p)-\ \a[{l-p)/p] on the NL, we easily recover (26). 

We can summarize the difference between lnC(p) and S w as follows. In lnC(p) we count 
(the logarithm of) the total number of different equivalence classes 7r that can be obtained for 
given p, i.e. that have P w > 0: lnC(p) thus measures the size of the support of the distribution 
P w . On the other hand, S n is the entropy of the distribution, and the fact that it is smaller 
than mC(p) indicates that the distribution is strongly (i.e. exponentially narrowly) peaked 
rather than uniformly spread over its support. 

4-5 Lower bound on the ground-state energy 

The inequality (25) implies the relation 

nH(p) > nH(p ) - N. (30) 

If we regard p as a function of po through the definition E g (po) = — n(l — 2p), then the above 
inequality gives a lower bound on p(po), or equivalently a lower bound on E g (p ). The result is 
shown in Fig. 7 together with numerical estimates of E g . Our lower bound is not particularly 
tight numerically, e.g. E g {p)/N > —1.56 at p = \ whereas numerically E g {p)/N is around 
— 1.40. However, this result has non-trivial significance because it is, as far as we know, the 
first analytical lower bound on the ground-state energy of the two-dimensional ±J Ising model 
with general p. 
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5. Summary and discussion 

We have derived upper and lower bounds on the number of equivalence classes of error 
patterns in the toric code. It has been shown that this number saturates its upper bound at 
an intermediate value of the error probability where we expect no singularities in physical 
quantities. This apparently non-conventional behaviour has been explained by noting that 
the number of typically realized equivalence classes, which is relevant for physical quantities, 
is significantly smaller than the number of equivalence classes with uniform weights given 
to all the cases. The logarithm of the former number has been shown to be equal to the 
information-theoretical entropy of the probability of error classes, which is further related to 
the thermodynamic entropy on the NL and therefore has a singularity (albeit a weak one) at 
the multicritical point. 

One of the upper bounds was compared with a lower bound, the latter involving the 
ground-state energy of the ±J Ising model, leading to a lower bound on the ground-state 
energy of this spin glass model. Although the resulting value of the lower bound is not nec- 
essarily impressive numerically, it is interesting that bounds on the number of equivalence 
classes of the toric code lead in a simple manner to a bound on the ground-state energy of a 
spin glass. The correspondence between the two problems was proposed by DKLP, 1 and we 
have exploited it here to derive an explicit result on a physical quantity. 

The present work would serve as a starting point for further developments based on the 
correspondence of two completely different problems. For example, an improved upper bound 
for C(p) will lead to a better lower bound on the ground-state energy. Improvements of the 
lower bound of C(p) may also be possible, but one should remember that such a result may 
not lead to an improved bound on the ground-state energy unless the obtained lower bound 
of C(p) is related to the latter quantity. 
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